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We study cosmological models derived from higher-order Gauss-Bonnet gravity F{R, G) by using 
the Lagrange multipher approach without assuming the presence of additional fields with the excep- 
tion of standard perfect fluid matter. The presence of Lagrange multipliers reduces the number of 
allowed solutions. We need to introduce compatibility conditions of the FRW equations, which im- 
pose strict restrictions on the metric or require the introduction of additional exotic matter. Several 
classes of F{R, G) models are generated and discussed. 

PACS numbers: 04.50.Kd, 04.20.Jb, 04.20.Cv, 98.80.Jk 



I. INTRODUCTION 

Astrophysical data indicate that the observed universe is in an accelerated phase |1|. This acceleration could be 
induced by the so-called dark energy (see Ref. |3| for a recent review and references therein) the nature and properties 
of which are not yet understood at fundamental level. 

There are several models for dark energy, but, in general, such a constituent can be figured out, in a coarse grain 
approach, as being constituted by some ideal fluid whose equation of state (EoS) exhibits non-standard properties, in 
particular, a non standard adiabatic index. On the other hand, dark energy can be considered as a global phenomenon 
associated with modifications of gravity [3], which are given by extensions of General Relativity (implying modifications 
of the Hilbert-Einstein action by introducing scalar fields, curvature invariants like the Ricci scalar R or the Ricci and 
Riemann tensors R^^, and -Ra^i^ or the Gauss-Bonnet topological invariant 0, Q)- Such extended actions describe 
effective theories coming from fundamental interactions [6| and, from a cosmological point of view, they should 
consistently describe the early-time inflation and late-time acceleration, without the introduction of any other material 
dark component. In this picture, the dark side of the universe could be traced as the lack of a final theory of gravity 
where gravitational interaction does not behave in the same way at all scales [7|. 

In fact. General Relativity presents problems at ultra-violet (Quantum Gravity) and infra-red (Cosmology) scales 
[a, |8| . A similar situation appears for dark matter phenomena. No fundamental candidate has been revealed up to 
now and dynamics of self- gravitating structures could be addressed by modifications of gravity like /(i?)-gravity (see 
for example [9]). A further problem is that it is not clear why dark energy had no effect at early epochs while it gives 
dominant contributions in today observed universe. According to the latest observational data, dark energy currently 
accounts for about 68 -^ 70 % of the total mass-energy amount of the universe (see, for example, the very recent data 
coming from the PLANCK collaboration [l^)- The main feature of dark energy is that its EoS parameter wd is 
negative: 



WD ^Pd/pd < 0, 



(1) 



where pu is the dark energy density and po the pressure. Due to this feature, the Hubble fluid is accelerating 
instead of decelerating as it should be if governed by ordinary perfect fluid matter. Although current data favor the 
standard ACDM cosmology, the uncertainties in the determination of the EoS dark energy parameter w are still too 
large, namely w — — 1.04lg'^g. Hence, one is not able to determine, without doubt, which of the three cases: w < —I, 
«; = — l,orw>— lis the one actually realized in our universe (lll.ll2|. Furthermore, the evolution of the cosmological 
constant is a very likely feature in order to connect the early quantum gravity regime with today observed universe 



Recently, a new cosmological model describing the different stages of evolution of the universe has been proposed 
[13 |. This model, named Dust of Dark energy, describes the accelerated expansion of the Universe by two scalar fields 
where one of them is a Lagrange multiplier (LM) which imposes a constraint on the dynamics [l4| - [la | . The exteirsioir 



to f{R) gravity via the addition of a Lagrange multiplier constraint has been proposed in Ref. [1J|. Such a model can 
be considered as a new version of modified gravity because dynamics and cosmological solutions are different from 
the standard version of f{R) gravity without such a constraint. This result is clear from a dynamical viewpoint: LMs 
are anholonomic constraints capable of reducing dynamics [17, 18]. Furthermore, using the LM approach helps in the 
formulation of covariant renormalizable gravity [19i] and can be re-conducted to the existence of Noether symmetries 
into dynamics P, [l6| ■ 

In the present paper, we study the Gauss-Bonnet gravity with LM constraints in view to recover realistic dark 
energy behaviors. Here we do not specify the form of the gravitational action but generically require a function of 
the Ricci scalar R and the Gauss-Bonnet invariant G, by considering a generic F{R, G')-gravity. In this way, due to 
the form of the Gauss-Bonnet invariant, all the curvature invariants are taken into account without arbitrary choices. 
Unlike the models dealt in [1^1, where the introduction of LM helps to generate a number of new solutions, we found 
the opposite effect in the case under consideration. Namely, we see that the appearance of LMs reduces the number of 
solutions and does not help in their generation. This feature is strictly related to the indetermination in the F{R, G) 
function endowed with all the gravitational degrees of freedom apart from the standard one in the Hilbert-Einstein 
action R. The paper is organized as follows. In Sect. II, we discuss the main features of F{R, G)-gravity where a LM 
is considered for the Gauss-Bonnet component. In view of studying dark energy models, the cosmological equations 
are derived. In particular, we recast the further gravitational degrees of freedom in the form of an effective fluid in 
order to derive a suitable EoS. Sect. Ill is devoted to the search for cosmological power law solutions for some classes 
of the F{R, G) function. Results are discussed in Sect. IV where we draw also our conclusions. 

II. F{R, G) GRAVITY WITH GAUSS-BONNET LAGRANGE MULTIPLIER 

As discussed in J3|, Ij, [l6| , higher-order and non local curvature invariants come out in any effective theory derived 
from unified approaches like strings, braneworld, etc. Here, we will consider a generic F{R, G) function of the Ricci 
scalar R and the Gauss-Bonnet invariant G, where we ask for a Lagrange multiplier for the Gauss-Bonnet term in 
order to constrain cosmology. The starting action has the following form: 

g If{R, G) + a (]^V^GV^G + F2{G)\ | + 5„, . (2) 

Here \/2k^ is the gravitational coupling constant, F{R,G) and F2{G) are arbitrary functions, Sm is the standard 
perfect find matter action, A is a Lagrange multiplier whose variation yields the constraint equation 

iv^GV'^G + F2(G) = 0. (3) 

G is the Gauss-Bonnet invariant defined as 

G = R^.apR'"'"^ - AR^.R'"' + i?2 . (4) 

The gravitational field equations are 

-^Fg^, + ^V^GV.G + F'j^R^, - V^V.F',^ + g^^OF'^^ + 2FL'q [R^^^^R.^fi^ + 2R^^p,R^^- 

— 2i?api?"^ + RRfiv) — 4 {RfiaPu + RiivQafi + RafiQiiV — Rafig^i^ — Raiy9Pii+ (5) 

Here F'j^ = dF{R, G)/dR and FL'g = dF{R, G)/dG + \dF2{G)/dG - V^(AV^G). In the 4-dimensional space-time, 
we have the following expression 



2 (i?p" '^ RuaPf + 2R^j_ai3uR°' — 2RafiR"y + RRfj,u) — -^Gg^i, — 0. 



Hence Eq. ([S]) can be rewritten as 



~2^9iii' + 2'^ ^iGV yG -f F'^Rf^y — V ^V^F'j^ + y^^DF^ + ^FL'^Gg^i, — 4 {Rf^apy + Rtii'gap+ /g\ 

+Rapgiiu — Raiig(3u — Raug^ii + 2-^ {gap.gpv — g^iugap)) V^V^fLq — k t^^ . 

Let us now consider a Friedmann-Robertson- Walker (FRW) universe with the flat spatial metric 

ds^ = -dt^ + a{tf Y^ {dx'f . {7) 

4=1,2,3 



Then from Eq. ^ one obtains 
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FRW equations are specified as soon as one defines the effective density and pressure 
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[F'rR + G FL'g - F) - XG^ - QHF'r - 2AH^F'L'g 
- {F'rR + G FL'g - F) + AHF'r + 2F' r 
16H (h + hA FL'g + SH^FL'g } ■ 



F', iP^ + 2;? 



(8) 



(9) 



(10) 



Here poff = 3k ^H^ and pcs = — « ^ (2H + 3H^] , respectively. A general method for solving such equations does 

not exist. However, a few solutions for general F{R, G) theory have been found |20l - [24| . We can reduce the problem, 
if we assume that 



FL'a = F'(., XG^^2^pM, Pm = 0. 



(11) 



We obtain an equation for the Lagrange multiplier, which will lead us to a model where the material field is in the 
form of dust (p = 0) . In such a case the conservation equation for matter 



p'u + 3iJ(pM + Pm) = , 



(12) 



could result violated but the general Bianchi identities have to remain conserved. Let us consider in more detail Eqs. 
([9|)- (|10p . It is easy to remove from these equations the constraint derived from the Lagrangian multiplier. We obtain 
then the compatibility condition of the equations, which takes the form 



SHk^Pm + 3i?K>M + K PM + 288HX 



2H 2H 



H 



HH 



(13) 



(h^ (-6H^ + 96H^ + H{-1 + 96H'^ + 24H^H] 

'2H^ + 3iJ(-l + 16H'^ + 8H^H)) H + QH^hA - H^H 



= 0. 



This equation can be rewritten in terms of the Gauss-Bonnet invariant obtaining 

AG (g^ - 2g) = -4^2 {3H{pM + Pm) + Pm) ■ 



(14) 



If the matter satisfies the conservation law we need to impose the additional conditions G = const or G = ci 
21n(i + C2). 

A similar situation arises when we consider another Lagrangian of slightly different type, that is 



S^^ fd''xV^lFiR,G) + x(^V^RV^R + F2iR)] }+Sm. 



(15) 



where A is again a Lagrange multiplier and F2 is an arbitrary function of the scalar curvature R. 

In this case the compatibility condition of the equations of motion will have exactly the same form, but with the 
replacement of the invariant of the Gauss-Bonnet with the Ricci curvature scalar 



AT? (7?2 - 2i?) = -4^2 {3H{pM + Pm) + Pm) 



(16) 

Clearly Eqs. (fT4|) and p6D do not contradict the validity of the general Bianchi identities for the theory but point 
out that (Gauss-Bonnet) curvature terms, selected by the Lagrange multipliers, can act as source terms. With 
these considerations in mind, let us search for cosmological solutions for F{R, G) gravity according to the Lagrange 
multiplier method. 



III. COSMOLOGICAL SOLUTIONS FOR F{R, G) GRAVITY 

Our aim is now to search for cosmological solutions compatible with F(i?, G) gravity that, eventually, select the 
form of the function. Let us start with power law solutions where the scale factor behaves as 

a(t) = aoi''", (17) 

so that the Hubble rate is given by 

F=^. (18) 



The FRW equations assume the form 

= (F-2K2pM)i"-24(/io-l)/io^t'"i^(i?,G)G-6(/io-l)M''i^(i?,G)i? + 

+ 6/io (-384(/io - l)/io^ (48(/io - l)/io^(9 + h^) + (29 - /io(14 + ih^))t^) A+ (19) 

+ [Ahfi't^°F{R,G)'(j + t^^F{R,G)'„ + ?,M{hf, - l)/io^ ((48(/io - l)/io^ + (9 - 2hQ)t^) A - i^X))) , 
and 
= (F + 2K2pM)t"-24(;io-l)/io^i'°i^(^,G)^ + 2(l-3/io)M''i^(-R,G)i^ + 

+ 2 (-1152(/io - l)/io^ (l6(/io - l)/io^(10 + /io)(3/io - 13) ~ (/lo - 3)(6 + /io)(3/io - 5)^") A+ (20) 

+ t (8{ho - l)ho^t^°F{R, G)'(. + 2hot^^F{R, G)'^ - 405504/io*^A + 847872/io^A - 479232/io^°A + 36864/io" A+ 
+ 42240/io^t^A - 64896Ct^A + 23808/10^^"^ - 1152/io^t*A + 4/io2t"F(i?''G)'g + t^^F{R]G)'p. + 18432/io^iA - 
- 36864/io^iA + 18432/io^°tA - 6528/io^i^A + 8448/io^t^A - 1920/io^i^A - 384(/io - l)/io^i^A(^)^ 

If we choose the Lagrange multiplier as 

A = ce ?5 ^5-3/10^ (21) 

were c is a constant, matter behaves as "dust" assuming the form 

9216c 12(ho-l)ho3 ,x9, K t; q^ „ ,„„^ 

PM^^T^e ^ (/lo - l)'Ci"'""''\ PM=0, (22) 

where the pressure is clearly zero. 

Choosing, instead, a Lagrange multiplier of the form 

A = Ait^^ , (23) 



one can use the reconstruction method proposed in Refs. |l6l l20l - |22| . However, in this case we need to introduce the 
matter in the special form 

, ^ , ,-3ho(i+«,) , 4608(feo - lfhl\,t-^^+>^- ( 48(/io - l)feo^ 5^^ A 

^^ ' «2 Ul4 + A2 + 3/io(l + ii;) -10 + A2 + 3/io(1 + H;' ^ ^ 

PM = wpM- (25) 

If A2 = 14 — 3/io(l + w) then the matter density takes the form 

^-3ho(i+») r + 1152(^0 - ifhlX, {-5t^ + 192(/io - l)feg Int) \ ^ ^^^^ 

On the other hand, if A2 = 10 — 3/io(l + w), the matter density takes the form 

/ 4608(feo - l)'/^gAi (_ i^(^'°, D^'o „ 5in^) \ 
PM = t-^''°(^+-) ci + ^ ^^A^ '- L . (27) 

With proper functions P{<j>), Z{<j>) and Q{4>) of a scalar field (p, which we can identify with the time t, we represent 
the term T{R, G) in the action in Eq. ([2]) as P{t)R + Z{t)G + Q{t). Varying this action with respect to t, we find 
{dP{t)/dt) R + {dZ{t)/dt) G + dQ{t)/dt = 0. By solving this equation, we obtain t = t{R, G). Combining this and 
using the above representation, we get the special case J^{R, G) = P{t)R + Z{t)G + Q{t). 

Let us consider now a few cases that can be worked out in detail. 



A. The case of F{G) gravity with a Lagrange multiplier 

Let us start with a function the form F = R + F{G), where we consider a pure Gauss-Bonnet function as correction 
to Einstein gravity. Considering ci = (without ordinary matter) then it is easy to find an exphcit form for F{G). 
There are two cases 



1. We have 



^ 4(feo ~ l)hl{3ho - 2) (-10 + ho + 3hl) 
^ " (/lo- 2)2(5 + /io) (-1 + 2/10 + 3/1^)^2 ' 

(2 - 3/to)tS 

1152/i4 (-2 + 7ho - hi - T/ig + 3h^) ' 

23/431/4 /;j3 + ;j4) 1/4 

t ^ ^i7i ^^, (28) 



F{G) 



2^|G3/2^(/io_l)/i3(20 - 32/io - 3/ig + 9h^) 
3{ha - 2)2/io(5 + Hq) (-1 + 2/io + 3hl) 



where the Gauss-Bonnet function scales as a 3/2 power. The dark energy behavior means that the power ho in 
Eq. pT|) is larger than 1. 



2. Another case is 



, _ (84 - 20/10 - 143h0^ - 30/ii^ -f 9/i;^) i*^ 

'^ ^ 192(-3 + ho){-l + ha)h^{2 + 3ha) (2 - 7ho + 3/ig) ' 

18432{-l + ho)^hl{2-7hn + 3hl)' ^ ' 



v/(-l + /io)/io (-84 + 20/10 + 143/ig + 3Qhl ~ %hfj 



VG{-3 + ho)ho{2 + 3ho) (2 - 7ho + 3hl) 
If we eliminate the Lagrange multiplier we get the following form of the function F{G) 

F{G) = p5i(y2 3i/4 + llGi/4)Gi/i2 (30) 

and then we get another fractional power for F{G). It is worth noticing that this kind of fractional behavior 
for generic gravitational actions is usually determined by the presence of Noether symmetries for the dynamical 
system (see Q for details). This shows the strict relation between LM method and the Noether Symmetry 
Approach to reduce dynamical systems. 

B. The case of F{R,G) gravity in presence of ordinary matter 

Let us assume now a generic F(R, G) function with no Lagrange multiplier, but the presence of matter in the form 
of a standard perfect fluid like pm = poa-^^^+'"\ pm = w poa-^'^^+'"\ The FRW Eqs. ^^ take the following form 

- 2K2pot4 ^ a3i3/io(ag^fto)3^(^^4 _ 6/io(4gi/i2(-i - g^ + ho)tf' + ri(-l + ho- r2)t2+'-^)) = 0, (31) 



aoH'"'"{aot''"f'" {~Ft^ - 8gi{g2 - ho)hl{-3 + /s + 3ho)tf- + 2n{ho - r2)(-l + 3ho + r2)t^+'-') - 

-2K'^pot^w = 0. (32) 



Here, the functions F^ and Fq can be selected in the form F'j^ = riV^ and Fq = git^^ . Using the reconstruction 
method ^^^ which we described above (that is Z{t) = F^, P{t) = F'^ and Q{t) = F{R,G) - F'^^R - F'^G), we 
can find t as a function of R and G, and then F as a function of R and G. However, this will lead to an equation 
containing an irrational degree for i, which cannot be solved in general. As above, one can consider a few special 
cases, but just for fixed values of h^ and w. In fact, deriving the explicit form of i^ as a function of R and G is very 
difficult. We choose, for the sake of simplicity, the case when F(i?, G) = Fi{R) +F2{G). Then the solution is possible, 
for po ^ 0, if 

1. .g2=4-3/io(l + w) 7^2 + r2 , 

2. 2 + r2=4-3/io(l + w)^ff2, 

3. 52 = 2 + r2 = 4 - 3/io(l + w) . 

Let us consider the first case. Solving Eqs. ((3T|) - ((32l) we find 



Po = -~4.al+^'"gihl{-4 + 3ho{l + w){~l + ho{4 + 3w))/{K,\l + w)), (33) 

1 



r2 = ^{l + ho±^l + 10h„ + hlj. (34) 

It is evident that, in order that the condition pQ > he satisfied, we need to impose restrictions on gi and ho, namely: 
gi > and 1/4 < ho < 4/3 or gi > and ho < 1/4 or ho > 4/3. Finally the functions Fi{R) and F2(G) assume the 
form 

„i (9+h„±^l + 10ho+hl) „Ul + ho±^l + 10ho+hl) ,, , , ^^A{'^+ha±^Jl + Wha+hl) Jy\U-hoT^/'^ + Who+hl) 

Fi{R) = — -^— ^jfto^^/ iQ Ijj V ^_K^ ^ 

3-hoT Vl + lO/io + hi 

F2(G) = 51^-^ — ^7 ^' ^ G'^''«(i+™). (36) 

In the second case the functions F12 assume the forms 

p,j^. ^ ri2^-i'-o(i+-)3-i^°(i+-)(2feo - l)(feo(2feo - l))-lMi+^) ^3,^(,^,^^ 

1 + w 

F2(G) = -gi22+^r^3^^(/io - 1)^(''«"^'/io^^G-"t. (38) 

If Po = then the function F is the sum of the functions F2 (G) of the first case ([36|) and the functions Fi (i?) of the 
second case ([37| . 

For the third case we easily find the function F by summing up the functions 

p,j^. ^ ^^ 2^-i'^o(i+^)3-i'-o(i+-)(2feo-l)(feo(-l + 2/^o))-^""^^+-) ^3,„(,^,^^ 

1 + w 

F2(G) = 31^ ^^^^^ ^ Gi'^^^'+^l (40) 

As general remark, also in these cases the functions F are fractional power laws pointing out that their form could 
be easily related to the presence of a Noether symmetry for the dynamical system. 



IV. DISCUSSION AND CONCLUSIONS 

We have examined cosmological models derived from F{R,G) gravity in presence of LMs like in the action ([2]). 
No other additional field has been considered with the exception of ordinary perfect fluid matter. Despite of Refs. 
|15l.[la|. where with the help of LMs new solutions were generated, here, we have constraints like (|14l) and P^ that 
have to be taken into account restricting the possible forms of the function F{R,G). On the other hand we can 
introduce, in addition to standard matter (in the form of a perfect fluid satisfying the conservation equation) , some 
form of exotic matter, which should compensate, in Eqs. (|14|) and ()16p . the terms corresponding to LMs. For a power 
law behavior of the cosmological scale factor a{t) and the LM of the form (|21l) . we have to introduce matter in the 
form (122), for which our model reduces to the standard F{R, G) model without LMs. 

To generate other solutions one can choose the LM in the form ([23|) ; then the matter of the form ((27|) will provide 
the compatibility of FRW equations. All of these examples are obtained for power law behaviors of the scale factor. 
Another example is based on F{G) gravity. In this case it is easy to reconstruct the explicit form of the function F{R). 
For a LM in the form © this function has the form ([251) and ([2^)) : with the exclusion of the LM, the function F{G) 
has the form (|5n|. All of these solutions have been built in the absence of ordinary matter, but the addition of such 
matter does not change much the situation. As examples, we have built several solutions for F{R) + F{G) gravity 
with matter, but without LMs like ([50]) - (|^0)) . Adding LM, as noted above, leads to the appearance of exotic matter 
and significant complication in equations: this means a series of additional restrictions on the F{R^ G) functions. 
In conclusion, the introduction of LMs can lead to additional restrictions, which can be avoided in the models with 
further fields, such as scalar fields. In this case, the compatibility of the FRW equations is achieved by the right 
choice of a scalar field. Standard choices for it, for example, = Ini, allow to give rise to consistent sets of equations 
as discussed in [16] . Finally, we can say that the LM method is a sort of inverse scattering approach that allows to 
reconstruct a model, given the cosmological behavior that one wants to fit. This feature is particularly relevant with 
respect to the issue of matching models with observational data instead of imposing, a priori, arbitrary theories. 
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